Single-cell transcriptomics of the human parasite Schistosoma mansoni first intra-molluscan stage reveals tentative tegumental and stem-cell regulators.

Schistosomiasis is a major Neglected Tropical Disease, caused by the infection with blood flukes in the genus Schistosoma. To complete the life cycle, the parasite undergoes asexual and sexual reproduction within an intermediate snail host and a definitive mammalian host, respectively. The intra-molluscan phase provides a critical amplification step that ensures a successful transmission. However, the cellular and molecular mechanisms underlying the development of the intra-molluscan stages remain poorly understood. Here, single cell suspensions from S. mansoni mother sporocysts were produced and sequenced using the droplet-based 10X Genomics Chromium platform. Six cell clusters comprising two tegument, muscle, neuron, parenchyma and stem/germinal cell clusters were identified and validated by in situ hybridisation. Gene Ontology term analysis predicted key biological processes for each of the clusters, including three stem/germinal sub-clusters. Furthermore, putative transcription factors predicted for stem/germinal and tegument clusters may play key roles during parasite development and interaction with the intermediate host.

Genome-wide transcriptome analysis of Echinococcus multilocularis larvae and germinative cell cultures reveals genes involved in parasite stem cell function.

The lethal zoonosis alveolar echinococcosis is caused by tumour-like growth of the metacestode stage of the tapeworm Echinococcus multilocularis within host organs. We previously demonstrated that metacestode proliferation is exclusively driven by somatic stem cells (germinative cells), which are the only mitotically active parasite cells that give rise to all differentiated cell types. The Echinococcus gene repertoire required for germinative cell maintenance and differentiation has not been characterised so far. We herein carried out Illumina sequencing on cDNA from Echinococcus metacestode vesicles, from metacestode tissue depleted of germinative cells, and from Echinococcus primary cell cultures. We identified a set of ~1,180 genes associated with germinative cells, which contained numerous known stem cell markers alongside genes involved in replication, cell cycle regulation, mitosis, meiosis, epigenetic modification, and nucleotide metabolism. Interestingly, we also identified 44 stem cell associated transcription factors that are likely involved in regulating germinative cell differentiation and/or pluripotency. By in situ hybridization and pulse-chase experiments, we also found a new general Echinococcus stem cell marker, EmCIP2Ah, and we provide evidence implying the presence of a slow cycling stem cell sub-population expressing the extracellular matrix factor Emkal1. RNA-Seq analyses on primary cell cultures revealed that metacestode-derived Echinococcus stem cells display an expanded differentiation capability and do not only form differentiated cell types of the metacestode, but also cells expressing genes specific for protoscoleces, adult worms, and oncospheres, including an ortholog of the schistosome praziquantel target, EmTRPMPZQ. Finally, we show that primary cell cultures contain a cell population expressing an ortholog of the tumour necrosis factor α receptor family and that mammalian TNFα accelerates the development of metacestode vesicles from germinative cells. Taken together, our analyses provide a robust and comprehensive characterization of the Echinococcus germinative cell transcriptome, demonstrate expanded differentiation capability of metacestode derived stem cells, and underscore the potential of primary germinative cell cultures to investigate developmental processes of the parasite. These data are relevant for studies into the role of Echinococcus stem cells in parasite development and will facilitate the design of anti-parasitic drugs that specifically act on the parasite germinative cell compartment.

A single locus determines praziquantel response in Schistosoma mansoni.

We previously performed a genome-wide association study (GWAS) to identify the genetic basis of praziquantel (PZQ) response in schistosomes, identifying two quantitative trait loci situated on chromosomes 2 and 3. We reanalyzed this GWAS using the latest (version 10) genome assembly showing that a single locus on chromosome 3, rather than two independent loci, determines drug response. These results reveal that PZQ response is monogenic and demonstrates the importance of high-quality genomic information.

A single locus determines praziquantel response in Schistosoma mansoni.

We previously performed a genome-wide association study (GWAS) to identify the genetic basis of praziquantel (PZQ) response in schistosomes, identifying two quantitative trait loci (QTL) situated on chromosome 2 and chromosome 3. We reanalyzed this GWAS using the latest (v10) genome assembly showing that a single locus on chromosome 3, rather than two independent loci, determines drug response. These results reveal that praziquantel response is monogenic and demonstrates the importance of high-quality genomic information.

The stage- and sex-specific transcriptome of the human parasite Schistosoma mansoni.

The flatworm Schistosoma mansoni is an important but neglected pathogen that causes the disease schistosomiasis in millions of people worldwide. The parasite has a complex life cycle, undergoing sexual reproduction in a mammalian host and asexual replication in a snail host. Understanding the molecular mechanisms that the parasite uses to transition between hosts and develop into dimorphic reproductively competent adults may reveal new strategies for control. We present the first comprehensive transcriptomic analysis of S. mansoni, from eggs to sexually naïve worms. Focusing on eight life stages spanning free-living water-borne and parasitic stages from both intermediate and definitive hosts, we have generated deep RNA-seq data for five replicates per group for a total of 75 data sets. The data were produced using a single approach to increase the accuracy of stage-to-stage comparisons and made accessible via a user-friendly tool to visualise and explore gene expression ( https://lifecycle.schisto.xyz/ ). These data are valuable for understanding the biology and sex-specific development of schistosomes and the interpretation of complementary genomic and functional genetics studies.

Corrigendum: Dynamics of host immune response development during Schistosoma mansoni infection.

[This corrects the article DOI: 10.3389/fimmu.2022.906338.].

A drug repurposing screen for whipworms informed by comparative genomics.

Hundreds of millions of people worldwide are infected with the whipworm Trichuris trichiura. Novel treatments are urgently needed as current drugs, such as albendazole, have relatively low efficacy. We have investigated whether drugs approved for other human diseases could be repurposed as novel anti-whipworm drugs. In a previous comparative genomics analysis, we identified 409 drugs approved for human use that we predicted to target parasitic worm proteins. Here we tested these ex vivo by assessing motility of adult worms of Trichuris muris, the murine whipworm, an established model for human whipworm research. We identified 14 compounds with EC50 values of ≤50 µM against T. muris ex vivo, and selected nine for testing in vivo. However, the best worm burden reduction seen in mice was just 19%. The high number of ex vivo hits against T. muris shows that we were successful at predicting parasite proteins that could be targeted by approved drugs. In contrast, the low efficacy of these compounds in mice suggest challenges due to their chemical properties (e.g. lipophilicity, polarity, molecular weight) and pharmacokinetics (e.g. absorption, distribution, metabolism, and excretion) that may (i) promote absorption by the host gastrointestinal tract, thereby reducing availability to the worms embedded in the large intestine, and/or (ii) restrict drug uptake by the worms. This indicates that identifying structural analogues that have reduced absorption by the host, and increased uptake by worms, may be necessary for successful drug development against whipworms.

Using ChEMBL to Complement Schistosome Drug Discovery.

Schistosomiasis is one of the most important neglected tropical diseases. Until an effective vaccine is registered for use, the cornerstone of schistosomiasis control remains chemotherapy with praziquantel. The sustainability of this strategy is at substantial risk due to the possibility of praziquantel insensitive/resistant schistosomes developing. Considerable time and effort could be saved in the schistosome drug discovery pipeline if available functional genomics, bioinformatics, cheminformatics and phenotypic resources are systematically leveraged. Our approach, described here, outlines how schistosome-specific resources/methodologies, coupled to the open-access drug discovery database ChEMBL, can be cooperatively used to accelerate early-stage, schistosome drug discovery efforts. Our process identified seven compounds (fimepinostat, trichostatin A, NVP-BEP800, luminespib, epoxomicin, CGP60474 and staurosporine) with ex vivo anti-schistosomula potencies in the sub-micromolar range. Three of those compounds (epoxomicin, CGP60474 and staurosporine) also demonstrated potent and fast-acting ex vivo effects on adult schistosomes and completely inhibited egg production. ChEMBL toxicity data were also leveraged to provide further support for progressing CGP60474 (as well as luminespib and TAE684) as a novel anti-schistosomal compound. As very few compounds are currently at the advanced stages of the anti-schistosomal pipeline, our approaches highlight a strategy by which new chemical matter can be identified and quickly progressed through preclinical development.

A manually curated annotation characterises genomic features of P. falciparum lncRNAs.

BACKGROUND: Important regulation occurs at the level of transcription in Plasmodium falciparum and growing evidence suggests that these apicomplexan parasites have complex regulatory networks. Recent studies implicate long noncoding RNAs (lncRNAs) as transcriptional regulators in P. falciparum. However, due to limited research and the lack of necessary experimental tools, our understanding of their role in the malaria-causing parasite remains largely unelucidated. In this work, we address one of these limitations, the lack of an updated and improved lncRNA annotation in P. falciparum. RESULTS: We generated long-read RNA sequencing data and integrated information extracted and curated from multiple sources to manually annotate lncRNAs. We identified 1119 novel lncRNAs and validated and refined 1250 existing annotations. Utilising the collated datasets, we generated evidence-based ranking scores for each annotation and characterised the distinct genomic contexts and features of P. falciparum lncRNAs. Certain features indicated subsets with potential biological significance such as 25 lncRNAs containing multiple introns, 335 lncRNAs lacking mutations in piggyBac mutagenic studies and lncRNAs associated with specific biologic processes including two new types of lncRNAs found proximal to var genes. CONCLUSIONS: The insights and the annotation presented in this study will serve as valuable tools for researchers seeking to understand the role of lncRNAs in parasite biology through both bioinformatics and experimental approaches.

The gut microbial metabolic capacity of microbiome-humanized vs. wild type rodents reveals a likely dual role of intestinal bacteria in hepato-intestinal schistosomiasis.

Increasing evidence shows that the host gut microbiota might be involved in the immunological cascade that culminates with the formation of tissue granulomas underlying the pathophysiology of hepato-intestinal schistosomiasis. In this study, we investigated the impact of Schistosoma mansoni infection on the gut microbial composition and functional potential of both wild type and microbiome-humanized mice. In spite of substantial differences in microbiome composition at baseline, selected pathways were consistently affected by parasite infection. The gut microbiomes of infected mice of both lines displayed, amongst other features, enhanced capacity for tryptophan and butyrate production, which might be linked to the activation of mechanisms aimed to prevent excessive injuries caused by migrating parasite eggs. Complementing data from previous studies, our findings suggest that the host gut microbiome might play a dual role in the pathophysiology of schistosomiasis, where intestinal bacteria may contribute to egg-associated pathology while, in turn, protect the host from uncontrolled tissue damage.

Genomic landscape of drug response reveals mediators of anthelmintic resistance.

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities.

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Dynamics of Host Immune Response Development During Schistosoma mansoni Infection.

Schistosomiasis is a disease of global significance, with severity and pathology directly related to how the host responds to infection. The immunological narrative of schistosomiasis has been constructed through decades of study, with researchers often focussing on isolated time points, cell types and tissue sites of interest. However, the field currently lacks a comprehensive and up-to-date understanding of the immune trajectory of schistosomiasis over infection and across multiple tissue sites. We have defined schistosome-elicited immune responses at several distinct stages of the parasite lifecycle, in three tissue sites affected by infection: the liver, spleen, and mesenteric lymph nodes. Additionally, by performing RNA-seq on the livers of schistosome infected mice, we have generated novel transcriptomic insight into the development of schistosome-associated liver pathology and fibrosis across the breadth of infection. Through depletion of CD11c+ cells during peak stages of schistosome-driven inflammation, we have revealed a critical role for CD11c+ cells in the co-ordination and regulation of Th2 inflammation during infection. Our data provide an updated and high-resolution account of how host immune responses evolve over the course of murine schistosomiasis, underscoring the significance of CD11c+ cells in dictating host immunopathology against this important helminth infection.

Genomic evidence of contemporary hybridization between Schistosoma species.

Hybridization between different species of parasites is increasingly being recognised as a major public and veterinary health concern at the interface of infectious diseases biology, evolution, epidemiology and ultimately control. Recent research has revealed that viable hybrids and introgressed lineages between Schistosoma spp. are prevalent across Africa and beyond, including those with zoonotic potential. However, it remains unclear whether these hybrid lineages represent recent hybridization events, suggesting hybridization is ongoing, and/or whether they represent introgressed lineages derived from ancient hybridization events. In human schistosomiasis, investigation is hampered by the inaccessibility of adult-stage worms due to their intravascular location, an issue which can be circumvented by post-mortem of livestock at abattoirs for Schistosoma spp. of known zoonotic potential. To characterise the composition of naturally-occurring schistosome hybrids, we performed whole-genome sequencing of 21 natural livestock infective schistosome isolates. To facilitate this, we also assembled a de novo chromosomal-scale draft assembly of Schistosoma curassoni. Genomic analyses identified isolates of S. bovis, S. curassoni and hybrids between the two species, all of which were early generation hybrids with multiple generations found within the same host. These results show that hybridization is an ongoing process within natural populations with the potential to further challenge elimination efforts against schistosomiasis.

Identification of anti-schistosomal, anthelmintic and anti-parasitic compounds curated and text-mined from the scientific literature.

More than a billion people are infected with parasitic worms, including nematodes, such as hookworms, and flatworms, such as blood flukes. Few drugs are available to treat worm infections, but high-throughput screening approaches hold promise to identify novel drug candidates. One problem for researchers who find an interesting 'hit' from a high-throughput screen is to identify whether that compound, or a similar compound has previously been published as having anthelmintic or anti-parasitic activity. Here, we present (i) data sets of 2,828 anthelmintic compounds, and 1,269 specific anti-schistosomal compounds, manually curated from scientific papers and books, and (ii) a data set of 24,335 potential anthelmintic and anti-parasitic compounds identified by text-mining PubMed abstracts. We provide their structures in simplified molecular-input line-entry system (SMILES) format so that researchers can easily compare 'hits' from their screens to these anthelmintic compounds and anti-parasitic compounds and find previous literature on them to support/halt their progression in drug discovery pipelines.

Population genomics of ancient and modern Trichuris trichiura.

The neglected tropical disease trichuriasis is caused by the whipworm Trichuris trichiura, a soil-transmitted helminth that has infected humans for millennia. Today, T. trichiura infects as many as 500 million people, predominantly in communities with poor sanitary infrastructure enabling sustained faecal-oral transmission. Using whole-genome sequencing of geographically distributed worms collected from human and other primate hosts, together with ancient samples preserved in archaeologically-defined latrines and deposits dated up to one thousand years old, we present the first population genomics study of T. trichiura. We describe the continent-scale genetic structure between whipworms infecting humans and baboons relative to those infecting other primates. Admixture and population demographic analyses support a stepwise distribution of genetic variation that is highest in Uganda, consistent with an African origin and subsequent translocation with human migration. Finally, genome-wide analyses between human samples and between human and non-human primate samples reveal local regions of genetic differentiation between geographically distinct populations. These data provide insight into zoonotic reservoirs of human-infective T. trichiura and will support future efforts toward the implementation of genomic epidemiology of this globally important helminth.

Correction: MiR-277/4989 regulate transcriptional landscape during juvenile to adult transition in the parasitic helminth Schistosoma mansoni.

[This corrects the article DOI: 10.1371/journal.pntd.0005559.].

Characterising genome architectures using genome decomposition analysis.

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high-quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.

Defining the early stages of intestinal colonisation by whipworms.

Whipworms are large metazoan parasites that inhabit multi-intracellular epithelial tunnels in the large intestine of their hosts, causing chronic disease in humans and other mammals. How first-stage larvae invade host epithelia and establish infection remains unclear. Here we investigate early infection events using both Trichuris muris infections of mice and murine caecaloids, the first in-vitro system for whipworm infection and organoid model for live helminths. We show that larvae degrade mucus layers to access epithelial cells. In early syncytial tunnels, larvae are completely intracellular, woven through multiple live dividing cells. Using single-cell RNA sequencing of infected mouse caecum, we reveal that progression of infection results in cell damage and an expansion of enterocytes expressing of Isg15, potentially instigating the host immune response to the whipworm and tissue repair. Our results unravel intestinal epithelium invasion by whipworms and reveal specific host-parasite interactions that allow the whipworm to establish its multi-intracellular niche.